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Abstract 

We construct a Hamiltonian for the generation of arbitrary pure states of the 

quantized electromagnetic field. The proposition is based upon the fact that a 

unitary transformation for the generation of number states has been already 

found. The general unitary transformation here obtained, would allow the 

use of nonlinear interactions for the production of pure states. We discuss 

the applicability of this method by giving examples of generation of simple 

superposition states. We also compare our Hamiltonian with the one resulting 

from the interaction of trapped ions with two laser fields. 
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The generation of pure nonclassical states of the quantized electromagnetic field is a 
central topic in quantum optics. Several schemes have been already proposed, based either on 
the non-unitary collapse of the state vector (through atom measurements) within micromaser 
environments 0,0] or by using a cavity QED unitary time-dependent interaction [Q]. Both 
approaches involve individual atoms interacting with a single-mode cavity field, what would 
demand extraordinary control in a generation experiment. It is therefore interesting to seek 
for alternative methods for the generation of nonclassical light. A significant advance was 
the proposal, for the first time, of a unitary generator for the number state \n) by Kilin and 
Horoshko 0]. They derived an expression for a Hamiltonian Hn such that Hn\0) = \n) and 
Hn\n) = |0), i.e., they obtained a Hamiltonian that generates the state \n) from the vacuum 
state |0), but only during particular interaction times [t^ = {m + l/2)7r]. This could be 
experimentally achieved by pumping conveniently prepared nonlinear media, as discussed in 
i- 

In this paper we propose a generalization of Kilin-Horoshko's procedure for the generation 
of an arbitrary pure state which for instance, can be expressed as a linear superposition 
of number states 

M 



m = T.Cn\n), (1) 

n=0 

where M is the maximum number of photons in the radiation field. For certain interaction 
times, we assume the existence of a Hamiltonian ^|^) such that |0) = |\E') and |\E') = 
|0). The Hamiltonian may be written in the following general form 

M M 

H\^) = E Cn (|0)(n| + |n)(0|) - Co E \k){k\ + PFP, (2) 

n=0 k=0 

where P = I ~ Y.k=o \k){k\ and F is any Hermitian operator. The sums of projectors 
Y^k=Q \k)(k\ are needed in order to compensate all the terms present in the superposition 
state in Eq.(0). For the sake of simplicity we have considered the coefficients Cn real. Our 
task now is to find an operator F such that we are able to delete the higher order powers of a 
and a) in order to obtain the simplest possible expression for We start by decomposing 
F in a convenient way 
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Now if we substitute Eq.(|^) and the expression for P into Eq.(0), after some algebra we find 
the conditions that the functions fi{k) = {k\fi{d'^d)\k) must obey in order to cancel most of 
the terms in i^i^). These conditions are 

/o(0) = Co; /o(fc) = -Co; /^(O) = C^; /™(A;) = 0, (4) 

for k = 1,2, . . . M and m = 1,2, . . . M. Those are necessary but not sufficient conditions for 
the determination of the functions /^(ata), what introduces a certain degree of arbitrarieness 
in their choice. A simple form for /o and fm involving a finite expansion in the annihilation 
and creation operators d and is 

/o(dtd) = Co 2 (l — 0.^0,^ J^{b)d) — 1 ; fm{ci^a) = Cm J-'{a'd); m ^ 0, (5) 

where 

M 

J^{d^d) =Y,Aia^ay. (6) 

1=0 

The condition /m(0) = Cm demands that jF(d"''d)|0) = |0), which means that we must 
have Aq = 1. Moreover, we can use the condition fm{k) = [see Eq.(^)] to determine the 
remaining coefficients Ai. The successive application of the function jF(dtd) onto the M 
number states |1), |2) . . . \M) gives rise to the following M coupled linear equations for the 



coefficients Ai: 



1 + Ai + A2 + ... + Am = 



1 + 2Ai + 4^2 + . . . + = 



1 + MAi + MM2 + ... + M*Mm = 0. (7) 

The solution of this set of equations always exists, what completely determines the functions 
fm- Therefore there is a specific set of numerical coefficients Ai for every value of M, and 



which are of the form: Ai = ai/M\, A2 = a2/M\, . . . Am = Qm/MI, with < |a2| < • • • < 
|aAf-i| < |«Af|- In particular, um = (— 1)*^/M!. This resuh is convenient for us because 
the powers of a^a will be multiplied by increasingly smaller coefficients, i.e., the relative 
importance of the higher-order terms will be consistently diminished. 
The complete and already simplified Hamiltonian will then read 

H\^^ = /o (ata) + ^ fj^(a^a) a" + a^"^ J^ia^a)] . (8) 



We would like to remark that the function /o(a^a) is important for the establishment of the 
match between our Hamiltonian and the "physical" interaction Hamiltonian in a realistic 
experiment. Here, we have made a choice for /q [see Eq.(^] that allows the inclusion of the 
Kerr effect in the generation process, as we are going to show. 

In what follows we discuss applications of this method. It would be interesting to generate 
a state exhibiting nonclassical properties such as squeezing, antibunching and sub-Poissonian 
character, for instance. These effects occur simultaneously in a single state, namely the 
binomial state ||^, which admits an expansion as the one in Eq.(|I]), having a finite number 
of number state (real) coefficients different from zero 



Cn = B 
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n! (M-n)! 

Its photon number distribution P„ = {B^Y is a binomial distribution. The binomial states 
are caracterized by two parameters: p being the probability of emission of a single photon, 
and M the maximum number of photons in the field. They interpolate between a number 
state \M) (containing M photons), as p — 1, and a coherent state \a) (with real amplitude 
a = y/pM) as p — s> and M 00, i.e., they belong to a class of "intermediate states" . Their 
nonclassical properties, of course, are strongly dependent on the values of the "interpolation 
parameters" p and M Generalizations of binomial states also include the squeezed 
coherent states [0. 

We start by showing an explicit calculation of the Hamiltonian for the particular case of 
M = 1, i.e., the state being generated would be the superposition of the vacuum state with 



the one-photon state |0) = Co|0) + Ci|l). In this case we have that J-'{a^a) = 1 + Aia^d. 
From the requirement (1 + Aid''a)\l) = [fourth condition in Eq.(^], we obtain Ai = —1. 
The Hamihonian in Eq.(R|) then becomes 



%> = (1-P)'/' 



1 - Ad'^d + 2 (d^d 



'1 — d^d)d + dUl — d^d 



(10) 



On the other hand, if we pump with a classical field Ec = Ee~^^^ +E*e^^^ (polarized along 
the y direction) a nonlinear medium characterized by linear and nonlinear susceptibilities 
X*^^^ and x''^^ respectively, the coupling of the signal (polarized along the x direction), the 
pump and the output fields will be described by the following Hamiltonian [^]: 

Hnl = X^ia^a + x^^^W? + +X^^La^^d^ + X^^lya^a\E\' + xfUE\^\E\^ 

+ (xSa^i? + X%ya^'E^ + + x^^lya^\E\'E + H.c.) . (11) 

Both the pump and signal are travelling waves propagating along the z axis. The expres- 
sion for our generator [in Eq.(|lOD] has the same form as the Hamiltonian in Eg . ([TTD , what 
allows, at least in principle, to know how to choose and prepare a nonlinear medium in 
such a way that we end up with the Hamiltonian in Eq. ([T0|) , as we are going to show. 
The linear and nonlinear susceptibilities have to assume specific values in order to make 
possible the correspondence between both Hamiltonians. For instance, there is no term in 
Eq.(p!OD proportional to a''"^, which means that we must have x^xlyy — 0- By comparing terms 
proportional to d\ d'^'^d and to d^a, we obtain 

xS(^o)i? + xSLl^l'^ = -Xil.^ = P'^^^ (12) 

and 



X. 



^I{E,) + x%y\E? = -xil. = -2(1 - py/'. (13) 



These relations among the different susceptibilities are of the same form as the ones found in 
reference [^. However, in our case we are generating a state which is the one-photon state 
coherently superposed to the vacuum, fact which is embodied in the "tuning" parameter 



(s) 

p. The susceptibihties xljii of course a feature of the chosen crystal and the first-order 
susceptibihties xl]'' can be changed by the apphcation of a static field Eq. We can control 
the generation procedure by adjusting the values of the pump field E and the static field 
Eq in order to satisfy the relations (|12|) and (|T^). From Eq.(^2|) we see that the probability 
of having one photon in the field, p, is proportional to the pump field amplitude E, as 
one would expect. We note that the conditions in Eq.(P!^ and connect the quantum 
superposition principle, represented by the coherent superposition of the vacuum with the 
one-photon state, with "macroscopic features" such as nonlinear susceptibilities in a crystal. 

Our general scheme allows the generation, in principle, of virtually any pure state of the 
quantized field through some kind of nonlinear interaction. As a second example, we now 
construct the generating Hamiltonian for the superposition of the vacuum state |0) with 
the two-photon state |2), or \ip) = Co\0) + C2I2). In this case the solution of the system of 
equations in Eq.(^ gives us Ai = —3/2, A2 = 1/2, and the Hamiltonian will read 

1 C r 

Hi^^ = Co 1- 5(a^a) + 4(a^a)^ - l{a^af + —I J^{a^a)a'^ + a^^J^{a^a) ; 

J V2 L 

3 1 

J='{a^a) = 1- -a^a + -{a^ay. (14) 

We note that because the target state above does contain even photon numbers only, 
there are fewer terms in the generating Hamiltonian, what means that additional conditions 
must be imposed on the relevant susceptibilities in the corresponding nonlinear Hamiltonian. 
More specifically, in Eq.([T^) we have terms of the type a^"^ (^a'^cij , which means that a fourth- 
order nonlinear susceptibility (x^^^) should take part in the generation process. Recently 
there have been developments both theoretical P] and experimental in processes involving 
five-wave mixing in fluids. This kind of medium has several advantages, such as a flexible 
geometry, for instance. However, due to its intrinsic isotropy, processes involving modes 
in co-linear propagation are forbidden. In general there are more stringent conditions over 
fluid media than in crystals. This of course may favour our scheme in the sense that we 
have only a few terms present in our Hamiltonian. 

We would expect that the Hamiltonian in Eq.® for the generation of binomial states (in 



the limits of p — > and M oo) should be somehow equivalent to Glauber's displacement 
operator D{a) = exp(Q;a''" — a*a). In fact we have 



^i'|*)^|a>|0) = lim 
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i.e., the application of either the appropriate Hamiltonian or £)(a) onto the vacuum 
state leads to the same (coherent) state It is also worth verifying that the Hamiltonian 
in Eq.(|D (in the limit of p — 1) is equivalent to Kilin-Horoshko's Hamiltonian IQ for the 
generation on number states. In this case = 6n,M, and 
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|0) = |M). 



(17) 



The difference between Kilin-Horoshko's Hamiltonian and ours with respect to the number 
state case rests on the functions jF(d^a). They do not need to contain higher order powers of 
{a) a). In fact a simpler function JF = 1 — d)a/M is enough for the number state generation 
1^. However, we shall stress that in a real experiment, a nonlinear medium has in general 
all its nonlinear susceptibilities excited as it is pumped. This means that it is convenient to 
have the contributions of all powers of d)a (up to the Mth) in the interaction Hamiltonian, 
as it does happen in our approach. Our generalization, then, is consistent with previously 
known results. 

The effective implementation of this technique would constitute a challenging experi- 
mental problem, starting from the design of appropriate media. It would be interesting a 
comment on an alternative physical system, other than the electromagnetic field, that could 
be suitable for the accomodation of such a nonlinear Hamiltonian. Recently there has been 
a great deal of interest in the generation of nonclassical states of motion of a single ion 



confined into an electromagnetic trap. There have been carried out successful experiments 



for the production of number, coherent and squeezed vacuum states of the ion's motion [jlO 
On the other hand, Raman-type excitation of a trapped ion via two laser fields produces an 
intrinsically nonlinear Hamiltonian, which is conveniently expressed (one-dimensional case). 



in the following form ]TT|: 



Hion = ::hQh{n;v){tVaf + H.c.; where f,{n;r]) = e'"^ '^H iuiL L 



Here rj = 27rao/A is the Lamb-Dicke parameter, being ao(= 1/^2?™^) the size of the 
ground state of the harmonic potential of natural frequency z/i, i7 = 1] 1^22/2(^^^21 —^l) is the 
system's effective two-photon Rabi frequency, f2j are the single-photon Rabi frequencies, UJ21 
the atomic transition frequency, and cul is the laser frequency. The detuning between the 
lasers is A = kui {k = 1,2, . . .). This decomposition in Eq.(|l^) is similar to the one we used 
in our Hamiltonian [Eq.(P)], the difference being that in the ion's case the functions fk{n) 
are already fully specified. For instance, if we consider 17 -C 1, and A; = 1, we may write 

A(n;r/)^e-'^'/2(^l-^ata^, (19) 

and the interaction Hamiltonian will read 



Q i 1 — —d^a I a — Q*d^ i 1 — —d^d 
\ 2 \ 2 , 



(20) 



The Hamiltonian above has the same form as the one in Eq.(|ToD, for the generation of the one- 
photon state {p —>■ 1). We note that despite the similarity between both Hamiltonians, they 
are not totally equivalent, because for that we must have r] = \/2, which is in contradiction 
with our assumption 77 <C 1, apart form other differences. Therefore it is not so evident the 
correspondence between the quantum state engineering we are here proposing, and highly 
nonlinear systems such as driven trapped ions. Nevertheless, the approach adopted here, 
based on fixed interaction times, should be in our opinion further investigated in order to 
seek for a clearer connection between our tailor-made Hamiltonian and a scheme of excitation 
of ions by laser beams. 

8 



We have presented here the exphcit construction of a Hamiltonian which is a gen- 
erator of nonclassical states of the quantized electromagnetic field. This extends the class of 
methods being considered for this purpose, drawing together, as particular cases, previously 
known generators such as the number state as well as the coherent state ones. Having a 
generation scheme based on unitary transformations in hands, new possibilities arise for the 
engineering of nonclassical light under feasible controlled circumstances, e.g., through non- 
linear interactions. Moreover, this scheme has the advantage of involving travelling waves, 
what obviously minimizes the destructive effect of dissipation. 
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